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Abstract. The typical atmospheric boundary layer (ABL) 
over the southeast (SE) Pacific Ocean is featured with a 
strong temperature inversion and a sharp moisture gradient 
across the ABL top. The strong moisture and temperature 
gradients result in a sharp refractivity gradient that can be 
precisely detected by the Global Positioning System (GPS) 
radio occultation (RO) measurements. In this paper, the Con- 
stellation Observing System for Meteorology, Ionosphere & 
Climate (COSMIC) GPS RO soundings, radiosondes and the 
high-resolution ECMWF analysis over the SE Pacific are an- 
alyzed. COSMIC RO is able to detect a wide range of ABL 
height variations (1-2 km) as observed from the radioson- 
des. However, the ECMWF analysis systematically underes- 
timates the ABL heights. The sharp refractivity gradient at 
the ABL top frequently exceeds the critical refraction (e.g., 
— 157N-unitkm -1 ) and becomes the so-called ducting con- 
dition, which results in a systematic RO refractivity bias (or 
called A -bias) inside the ABL. Simulation study based on 
radiosonde profiles reveals the magnitudes of the A -biases 
are vertical resolution dependent. The A -bias is also the pri- 
mary cause of the systematically smaller refractivity gradient 
(rarely exceeding — 1 lON-unitkm -1 ) at the ABL top from 
RO measurement. However, the A -bias seems not affect the 
ABL height detection. Instead, the very large RO bending 
angle and the sharp refractivity gradient due to ducting al- 
low reliable detection of the ABL height from GPS RO. The 
seasonal mean climatology of ABL heights derived from a 
nine-month composite of COSMIC RO soundings over the 
SE Pacific reveals significant differences from the ECMWF 
analysis. Both show an increase of ABL height from the shal- 


low stratocumulus near the coast to a much higher trade wind 
inversion further off the coast. However, COSMIC RO shows 
an overall deeper ABL and reveals different locations of the 
minimum and maximum ABL heights as compared to the 
ECMWF analysis. At low latitudes, despite the decreasing 
number of COSMIC RO soundings and the lower percentage 
of soundings that penetrate into the lowest 500-m above the 
mean-sea-level, there are small sampling errors in the mean 
ABL height climatology. The difference of ABL height cli- 
matology between COSMIC RO and ECMWF analysis over 
SE Pacific is significant and requires further studies. 


1 Introduction 

The atmospheric boundary layer (ABL) over the southeast 

(SE) Pacific Ocean is characterized by a strong temperature 
inversion layer across the ABL top. Persistent shallow stra- 
tocumulus (Sc) cloud decks trapped below the ABL inversion 
play an important role in cloud-climate feedback processes 
because of its radiative cooling effect (e.g., Randall et al., 
1984; Slingo, 1990; Klein and Hartmann, 1993; Ma et al., 
1996; Bony and Dufresne, 2005; Clement et al., 2009). Lack 
of understanding and observations of complicated ABL pro- 
cesses (e.g., cloud radiative cooling, entrainment, turbulent 
mixing, ocean-atmosphere coupling) leads to great model 
uncertainties in weather and climate predictions (e.g., Ran- 
dall et al., 1998; Zeng et al., 2004; Palm et al., 2005; Jor- 
dan et al., 2010). Observations of thermodynamic struc- 
ture of the ABL over the open oceans are limited in space 
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and time to ground-based measurements over islands, ship 
soundings and several field campaigns. Conventional passive 
microwave or infrared satellite sounders do not have a suffi- 
cient vertical resolution to resolve ABL vertical structures, 
especially the thin transition layer across the ABL inversion 
layer. Several satellite measurements have been used to in- 
directly infer the ABL height through aerosol layer or cloud- 
top-height observations, which could serve as good tracers 
for ABL top over the well-mixed boundary layer regime 
(e.g., stratocumulus-topped ABL). The Cloud-Aerosol Lidar 
and Infrared Pathfinder Satellite Observations (CALIPSO) li- 
dar provides high precision height measurements of aerosol 
and low cloud layers (e.g., Wu et al., 2008). The Multi- 
angle Imaging SpectroRadiometer (MISR) on Terra satellite 
can retrieve the cloud-top-height (CTH) using a stereopho- 
togrammetric technique (e.g., Garay et al., 2008; Harshvard- 
han et al., 2009). Cloud-top-temperature (CTT) measure- 
ments from the geostationary or polar-orbiting satellites (e.g., 
Moderate Resolution Imaging Spectroradiometer, MODIS) 
through the infrared channel are also used for CTH re- 
trieval with supplemental information of the vertical temper- 
ature structures (Wood and Bretherton, 2004; Zuidema et al., 
2009). Nonetheless, all the above methods rely on the detec- 
tion of ABL tracers (e.g., aerosol or clouds) and encounter 
significant challenges over the broken cloud (e.g., trade cu- 
mulus) or clear-sky regions. Also, the requirement for a- 
priori meteorology information makes the MODIS CTH re- 
trieval somewhat model dependent. 

As a promising technique for profiling the shallow ABL, 
Constellation Observing System for Meteorology Ionosphere 
and Climate (COSMIC) radio occultation (RO) observations 
(~2000 daily profiles, available since April 2006) can mea- 
sure ABL properties with global coverage, high vertical 
resolution (^200-m in the lower troposphere) and cloud- 
penetrating capability. It is important to note that the imple- 
mentation of the open-loop tracking technique on the COS- 
MIC RO receivers significantly improves the quality of RO 
soundings in the moist ABL, as compared with the early 
GPS RO missions that used phase-locked loop (PLL) track- 
ing technique (Sokolovskiy 2001; Ao et al., 2009). Sev- 
eral studies have demonstrated the values of RO soundings 
in detecting the ABL height (e.g., Sokolovskiy et al., 2006; 
Sokolovskiy et al., 2007; Ao et al., 2008; Basha and Ratnam, 
2009; Guo et al., 2011; Ao et al., 2012). However, probing 
the ABL interior with RO remains to be challenging. One 
of the critical problems is the systematic negative refractiv- 
ity biases (A-bias) in RO sounding profiles, which is often 
pronounced over the subtropical eastern oceans (Xie et al., 
2010a), where the persistent sharp temperature and moisture 
gradient at the ABL top often results in the so-called ducting 
(or super-refraction) conditions (e.g., von Engeln and Teix- 
eira, 2004; Lopez, 2009). In particular, over the SE Pacific, 
the overall quality and the scientific value of the RO sound- 
ings inside the ABL are yet to be explored. 
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Fig. 1. Map of the ship-home radiosonde (circle) and COSMIC RO 
(cross) sounding locations during VOCALS-REx field campaign 
from 20 October to 1 December 2008. 

In this study, we analyzed the quality of COSMIC RO 
soundings in the SE Pacific by comparing the ABL top height 
measurements with those from the VAMOS Ocean-Cloud- 
Atmosphere-Land Study Regional Experiment (VOCALS- 
REx) radiosonde (Wood et al., 2011) and the high-resolution 
European Centre for Medium-Range Weather Forecasts 
(ECMWF) analysis data. The paper is organized as follows: 
Sect. 2 describes the datasets used in this study, followed by 
the methodology and validation results in Sect. 3; Sect. 4 
compares several key ABL parameters among COSMIC RO, 
VOCALS radiosonde observations and ECMWF analysis; 
Sect. 5 presents a seasonal mean ABL height climatology 
as well as the transect along ~20° S derived from COSMIC 
RO and ECMWF analysis; Sect. 6 discusses advances and 
limitations of current RO soundings for ABL study; Finally, 
conclusions and future work are given in Sect. 7. 

2 Data description 

This study focuses on GPS RO boundary layer measure- 
ments over the SE Pacific (Fig. 1) where the VOCALS-REx 
field campaign took place during 20 October-1 December 
2008 (Wood et al., 2011). The radiosondes collected from 
research cruises along the 20° S transect over SE Pacific 
in October 2006 and 2007 are also analysed (de Szoeke et 
al., 2010). The launch locations of the VOCALS radioson- 
des from NOAA research vessel (R/V Ronald H. Brown) 
and the COSMIC RO soundings during the campaign pe- 
riod are shown. A total of 213 radiosondes (10-m verti- 
cal sampling) were obtained at a 4-h interval. The COS- 
MIC RO soundings are provided by Jet Propulsion Labo- 
ratory (JPL), available at http://gnenesis.jpl.nasa.gov (e.g., 
Hajj et al., 2002; Ao et al., 2009). The spatial and tem- 
poral sampling of COSMIC RO soundings is described in 
Sect. 5. The refractivity retrieval is reported as a function 
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of geometric height above mean-sea-level (a.m.s.L), and the 
bending angle is reported as a function of impact parame- 
ter, which is the product of refractive index and the local 
curvature radius at the tangent point. The detailed descrip- 
tion of GPS RO occultation technique has been described 
in earlier studies (e.g., Kursinski et al., 1997; Anthes et 
al., 2008 etc.). Although RO soundings could potentially 
achieve ~60 m vertical resolution (Gobumov, 2004), the res- 
olution of JPL bending angle and refractivity profiles in the 
lower troposphere is limited by a ~200-m filter applied in 
the standard retrieval (similarly for the other RO data cen- 
ters, e.g., Ho et al., 2009). Since the marine ABL is the main 
focus of this paper, only RO profiles that penetrate below 
m.s.l. = 0.5 km are used. This screening eliminates almost 
all RO profiles over the continents. We also examine the 
high-resolution ECMWF global analysis data (TL799L91), 
which have a horizontal resolution of ~25 km (0.25° latitude 
x0.25° longitude) with 91 vertical layers from the surface to 
0.01 hPa (Bauer et al. 2006). The ECMWF six-hourly out- 
puts (00:00/06:00/12:00/18:00 Z) have 17 unevenly-spaced 
layers below 800 hPa (^2 km), with denser sampling near 
the surface. The vertical sample interval is ^200-m between 
900 and 800 hPa, decreasing to ^40-m near the surface. 

The ECMWF assimilation has incorporated the COS- 
MIC RO soundings in its analysis (Healy, 2008). However, 
ECMWF does not assimilate the RO profiles below the duct- 
ing layer at the ABL top (Poli et al., 2010), which is fre- 
quently seen in the analysis over subtropical eastern oceans, 
especially over SE Pacific (e.g., von Engeln and Teixeira, 
2004; Lopez, 2009). Also the VOCALS radiosondes are not 
assimilated into the analysis. Therefore, to large extent, the 
RO, the radiosonde and the analysis data can be considered 
as independent measurements in the ABL. 


3 Methodology and validation 

Two important ABL parameters are analyzed here, i.e., the 
ABL top height and the refractivity gradient near the ABL 
top. The ABL height is generally marked as a transition 
height with a strong temperature inversion and a sharp mois- 
ture gradient. It is a key parameter to characterize the ABL, 
which is determined by the vertical scale of turbulent eddies 
and controls flux exchanges among turbulent mixing, heat 
and moisture in the ABL. It also depends on dynamic con- 
ditions of the underlying ocean surface and the free tropo- 
sphere. The inversion strength, i.e., the temperature gradient 
across the thin transition layer (100-200 m thick) at the ABL 
top, is a measure of the stability of the transition layer. The 
inversion strength along with turbulent kinetic energy (TKE) 
controls the entrainment rate, a critical cloud dynamic pa- 
rameter that determines the magnitude of warming and dry- 
ing of the ABL by mixing its air with the free troposphere 
(Lilly 1968; Deardorff, 1976). The ABL inversion strength 
plays a critical role in evolution of cloud thickness and cov- 


erage (Wood and Bretherton, 2004). Although GPS RO can- 
not infer the ABL temperature inversion directly, the RO re- 
fractivity gradient depends strongly on both temperature and 
moisture gradients, and thus can be used as a valuable ob- 
servational constraint to the ABL inversion strength and the 
ABL top entrainment rate. 

3.1 Typical ABL structure over VOCALS region 

Figure 2a shows a typical ABL temperature and humidity 
structure from radiosonde observation over the VOCALS re- 
gion. The ABL top height is well defined by separating 
a moist, well-mixed boundary layer from a dry free tropo- 
sphere above. The stratocumulus cloud is trapped underneath 
the inversion base (at ~1.5 km) and is marked by the box 
with relative humidity (not shown) exceeding 94 % (Wang et 
al., 1999). A very interesting kink is seen in the radiosonde 
temperature profile at the cloud base (~ 1 km), showing a 
change in the temperature lapse rate from a nearly dry- 
adiabatic (~9.5 °Ckm -1 , a well-mixed ABL) to a saturated- 
adiabatic lapse rate (~4.3 °Ckm _1 ) inside the cloud. The 
near-coincident ECMWF analysis profile (Fig. 2b) shows a 
lower inversion base (i.e., cloud top height) at ~ 1.1 km with 
a smoother transition across the ABL top in temperature and 
specific humidity. The reduced sharpness is likely due to 
the coarser vertical resolution in ECMWF analysis (^200- 
m vertical sampling). As expected, the small structures of 
the transition from dry-adiabatic to saturated-adiabatic in the 
radiosonde are not present in the ECMWF profile. 

Based on the atmospheric temperature, humidity and pres- 
sure of the radiosonde and ECMWF data, the refractivity can 
be derived (Smith and Weintraub, 1953) and directly com- 
pared with COSMIC RO observations. It is important to 
note that the opposite changes in temperature and moisture 
across the ABL top lead to an enhanced decrease in refrac- 
tivity (Fig. 2c). Both the ECMWF and COSMIC refractivity 
profiles agree well with the radiosonde above the inversion 
base (~1.5 km). Below this level, the ECMWF profile shows 
a much reduced gradient across the ABL top, yet, matches 
well with the radiosonde below ~1.1 km. COSMIC refrac- 
tivity profile, on the other hand, shows systematically nega- 
tive errors inside the ABL. 

Assuming a local spherically symmetric (or horizontally 
homogeneous) atmosphere, one can compute the bending 
angle profile (Fig. 2d) as would be observed by GPS RO 
from a radiosonde or ECMWF refractivity profile (see de- 
tails in Sect. 4.3). The maximum bending angle generally 
corresponds to the MRG. As seen in Fig. 2c and d, the 
sharp refractivity gradient at the ABL top in the radiosonde 
(~1.5 km) is mapped into a sharp increase in bending angle 
at ~3 km (in impact height, i.e., the impact parameter mi- 
nus the Earth’s radius). The high-resolution COSMIC RO 
bending (5-m smoothed) captures the transition height with 
a large bending increase very well. The standard RO bend- 
ing (with the ^200-m vertical smoothing) also shows a large 
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Fig. 2. Typical ABL structure in one VOCALS radiosonde (black) with the near-coincident COSMIC-RO (blue) and ECMWF analysis 
(green) profiles, (a) Radiosonde specific humidity and temperature; (b) ECMWF analysis specific humidity and temperature; The boxes 
in (a, b) indicate the cloud region with relative humidity exceeding 94 %; (c) refractivity of radiosonde (black), COSMIC RO (blue) and 
ECMWF analysis (green) as well as simulated refractivity retrieval (black-dotted) based on the simulated radiosonde bending; the ABL 
heights are listed in the upper right comer, respectively; (d) simulated bending angle of radiosonde (black), ECMWF analysis (green) as well 
as the standard (blue-solid) and high-resolution (cyan) COSMIC RO bending angles. 


increase across the ABL top but reveals a smaller and lower 
peak bending, which results in a slightly lower MRG height 
(~ 1 .43 km). Although the ECMWF profile produces a much 
weaker ducting layer at a lower MRG height (~1.3 km), the 
impact height of the peak bending in this ECMWF profile is 
close to the radiosonde. Nevertheless, the bending below the 
peak is systematically larger than the radiosonde. 

Several studies have utilized such features in RO refrac- 
tivity profile for detecting the ABL height by locating the 
height of maximum refractivity gradient (MRG) in GPS RO 
profiles with some differences in algorithm implementation 
(e.g., Sokolovskiy et al., 2006; Ao et al., 2008; Basha and 
Ratnam, 2009; Guo et al., 2011; Ao et al., 2012). For a 
well-defined ABL (e.g., over VOCALS region), the refrac- 
tivity gradient method is consistent with other conventional 
definitions (e.g., gradient method based on temperature or 
humidity profiles) (Basha and Ratnam, 2009; Seidel et al., 
2010; von Engeln and Teixeira, 2011; Ao et al., 2012). How- 
ever, over the region without a well-defined ABL transition 
layer (e.g. the deep convection zone), the MRG method could 
result in significant discrepancies in ABL height comparing 
with other definitions (e.g., von Engeln and Teixeira, 2011). 
In this paper, we choose the MRG method and define the 


ABL height as the altitude where the refractivity gradient 
reaches its most negative value anywhere below 5 km. As the 
elevated ABL height is the main focus of this study, the low- 
est 50-ma.m.s.l. are excluded. This is because some strong 
moisture gradients found near the warm ocean surface in 
ECMWF analysis could significantly lower the ABL heights 
derived from the MRG method (von Engeln and Teixeira, 
2011). Note that the three datasets are associated with dif- 
ferent vertical resolution: radiosonde (10-m), RO (^200- 
m) and ECMWF (varying from~40-m to 200-m). To en- 
sure the same ABL height quality derived with the MRG 
method, we interpolate RO and ECMWF refractivity pro- 
files (cubic-spline) on a 10-m vertical grid before the gra- 
dient calculation. RO shows a reasonably good agreement 
with radiosonde in ABL height with a small negative error of 
^90-m, while ECMWF analysis has a larger negative error 
of ~200-m (Fig. 2c). 

Given the simulated radiosonde bending angle profile in 
Fig. 2d, the radiosonde refractivity profile expected from RO 
observation can be derived (dotted-line in Fig. 2c) following 
the standard Abel-inversion integration (e.g., Fjeldbo et al., 
1971; Kursinski et al., 1997). The simulated refractivity pro- 
file also shows negative errors against the original radiosonde 
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Fig. 3. (a) Scatter-plot of the ABL heights for near-coincident radiosonde (100-m smoothing) and COSMIC RO soundings with the linear fit 
(blue); (b) scatter-plot of the ABL heights for near-coincident radiosonde and ECMWF with the linear fit (green). 


profile. Such negative errors are solely caused by the so- 
call ducting phenomenon (i.e., refractivity gradient exceeds 
the critical refraction, — 157N-unitkm _1 ) at the top of the 
ABL where the Abel-inversion integration becomes a non- 
unique problem (Sokolovskiy, 2003; Ao et al., 2003; Xie, et 
al., 2006; Xie et al., 2010a). In this example, the COSMIC 
refractivity closely matches the simulated Abel-retrieved ra- 
diosonde refractivity profile, which strongly suggests the 
negative errors in COSMIC refractivity are primarily caused 
by the ducting. Additional negative errors in the RO refrac- 
tivity profile, relative to the radiosonde Abel profile, could be 
attributed to other factors such as the space and time differ- 
ence between the RO and the radiosonde (^268 km, ~1 h) 
and the representativeness difference between the point mea- 
surement for radiosonde and the limb sounding (horizontally 
averaged) observation for GPS RO. 

3.2 ABL height comparisons 

Out of a total of 190 radiosonde profiles from the VOCALS 
field campaign (15-25° S, 70-86° W), only 25 were found to 
be near-coincident with the quality-screened COSMIC RO 
profiles within 3h in time and 300 km in space. Further 
searching over the radiosondes (142 profiles) from cruises in 
October of 2006 and 2007, another 18 near-coincident pairs 
were found. Figure 3a compares the ABL heights for the 43 
matched cases (mean time difference ~1.5h; mean-distance 
^226 km), and a low bias (~ 100-m) is evident in the RO 
measurements. Linear regression analysis shows a high cor- 
relation coefficient (~0.85) with small 1 -sigma uncertainty 
(~90m). The limited number of near-coincident cases and 
the relatively large space/time difference between RO and 
the radiosondes could contribute to the ABL height differ- 
ence. The radiosonde-ECMWF (190 pairs) comparison, as 
shown in Fig. 3b, yields a low bias in the ECMWF analysis. 
The correlation is lower (M).54) and the 1 -sigma uncertainty 
is ~70 m. The systematic low bias in ECMWF ABL height is 
consistent with previous weather and climate model studies 


over the region (e.g., Bretherton et al., 2004; Rahn and Gar- 
reaud, 2010; Wyant et al., 2010; von Engeln and Teixeira, 
2011; Wang et al.,2011). 

4 Atmospheric boundary layer over southeast Pacific 
Ocean 

The short-period VOCALS campaign and the relatively 
sparse sampling of COSMIC RO soundings greatly limit the 
number of near-coincident pairs. COSMIC RO has a random 
sampling pattern especially in the zonal direction (see de- 
tails in Sect. 5), whereas radiosondes and ECMWF analysis 
have relatively uniform and dense sampling over the cam- 
paign area. An alternative approach is to evaluate ABL pa- 
rameters through their PDF (probability density function or 
histogram) using all available data sets over the region of in- 
terest. These PDFs can also yield differences among COS- 
MIC RO, radiosondes and ECMWF analysis. 

4.1 ABL height comparisons 

In this comparison the ABL heights are first derived from the 
refractivity profiles using the MRG method, and then sorted 
into 100-m bins between the surface and 4 km. The PDF is 
derived in each bin, normalized by the bin size and the total 
number of samples within 0-4 km height range. The PDF 
of the ABL heights for the 25 radiosondes (100-m smooth- 
ing) and the near-coincident COSMIC soundings as well as 
the ECMWF analysis are shown in Fig. 4a. The PDF of the 
COSMIC shows a better agreement with the radiosonde than 
the ECMWF. Both the RO and the radiosonde show a maxi- 
mum frequency of the ABL heights at ~1.5 km and exhibit a 
similar distribution. The ECMWF analysis, however, shows 
a lower ABL height at its peak frequency (~1.1 km), with 
the PDF shifting to the lower heights. 
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Fig. 4. The PDF of the ABL heights (a) for radiosondes (25 profiles, with 100-m smoothing) and the near-coincident COSMIC RO and 
ECMWF analysis; (b) for radiosonde (190 profiles) with three vertical smoothing intervals (100, 200 and 300 m) and the near-coincident 
ECMWF analysis over the region (15-25° S, 70-86° W) from 25 October to 1 December 2008; (c) for COSMIC RO soundings (70 profiles) 
and the near-coincident ECMWF analysis over the region (20-25° S, 70-86° W) from 25 October to 1 December 2008. For better illustration, 
the normalized PDF is divided by 3 to scale the maximum PDF close to 100. The shape of the PDF is not changed due to the scaling. 


A reduced vertical resolution results in a decrease in the 
vertical refractivity gradient and thus affects the ABL height 
derived from the MRG method. Figure 4b shows the PDF of 
the ABL heights for the 190 radiosondes (within 15-25° S) 
after applying three vertical smoothing levels (100, 200 and 
300 m). The ABL height at peak frequency (~1.3 km) does 
not change due to different levels of smoothing. The 200-m 
smoothing result is almost identical to the 100-m smooth- 
ing. While the 300-m smoothing does reduce the peak fre- 
quency of ABL height at 1.3 km, but the overall shape (e.g., 
the width) of the PDF is maintained. This is because the 
ABL top height is very well-defined, marked by one domi- 
nant sharp gradient in the lower troposphere, and is less sen- 
sitive to the smoothing. All three curves show the major- 
ity of ABL heights varying from 0.8 km to 2.1 km. On the 
other hand, the near-coincident ECMWF profiles (190 pairs, 
mean time difference = 1 .5 h; mean distance = 85 km) show a 
systematically lower ABL height (1.1 km) at the peak occur- 
rence frequency comparing to the radiosondes (Fig. 4b). This 
low bias (^200-m) in the ECMWF ABL height is consistent 
with the scatter plot in Fig. 3b. 

We also compute the PDF of ABL heights from COS- 
MIC soundings over the region (20-25° S, 70-86° W) dur- 
ing the VOCALS campaign period along with the near- 
coincident ECMWF analysis (70 pairs, mean time differ- 
ence = 1.5 h; mean distance = 130 km). COSMIC soundings 
show the peak frequency of ABL height at 1.3 km, while the 
ECMWF analysis shows systematically lower ABL heights 
with a peak occurrence frequency at 1.1 km (Fig. 4c). COS- 
MIC soundings are capable of resolving higher ABL with 
realistic range of ABL heights variation (0.9-2. 1 km) as re- 
vealed by the radiosondes (Fig. 4b). On the other hand, the 
ECMWF analysis systematically underestimates the ABL 
heights and shows limited capability of correctly reproduc- 
ing the ABL height variations. 


4.2 Maximum refractivity gradient (MRG) at the ABL 
top 

Sharp temperature inversions and drops in moisture both con- 
tribute to the large gradient of refractivity across the ABL 
top. Note that the inversion strength is an important ABL 
parameter that measures the stability of the transition layer 
and controls the moisture and momentum fluxes exchange 
between the boundary layer and the free troposphere. It also 
strongly affects the cloud-top entrainment rate and the evolu- 
tion of stratocumulus clouds below the inversion (Wood and 
Bretherton, 2004). The refractivity gradient measurement, 
therefore, is important for understanding the cloud-top en- 
trainment and other ABL processes. 

In this section we investigate the PDF of MRG at the ABL 
height using the near-coincident radiosondes, ECMWF anal- 
ysis and COSMIC RO as shown in Fig. 4b and c. To account 
for the large range of MRG values, we take the logarithm 
of the absolute value of MRG to base 10, and calculate the 
PDF with even bin size of 0.2 from 0 to 3. These bins are 
later converted back to the corresponding MRG from — 1 to 

— 1000 (N-unitkm _1 ) in linear scale. Again, the PDF that 
measures the number of occurrences in each bin is normal- 
ized by the bin size and the total number of samples. 

As expected, the refractivity gradient at the ABL top, or 
the MRG value, depends on the vertical resolution of the 
measurements. Lower vertical resolution leads to a smaller 
gradient. When the MRG is larger (more negative) than 
the critical refraction (— 157N-unitkm -1 ) the ducting con- 
dition occurs. With the original 10-m sampling, the MRG 
values of the radiosondes (190 profiles) vary from —200 to 

— 1000 (N-unitkm _1 ). To investigate the smoothing effect 
on the MRG, we apply different smoothing lengths on the 
radiosonde data. Figure 5a shows the PDF of MRG for the 
radiosondes and the near-coincident ECMWF. The PDF peak 
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Fig. 5. (a) PDF of MRG for the radiosondes (190 profiles, with different vertical smoothing) and the near-coincident ECMWF analysis over 
the region (15-25° S, 70-86° W); (b) PDF of the MRG for COSMIC RO soundings (70 profiles) and the near-coincident ECMWF analysis 
over the region (20-25° S, 70-86° W) from 25 October to 1 December 2008. The vertical black dashed lines denote the critical refraction 
gradient (— 157N-unitkm _1 ). The normalized PDF is multiplied by 20 for better illustration. 


shifts to a smaller MRG value when the smoothing length ap- 
plied on the radiosondes increases. After the 100-m or 200-m 
vertical smoothing, the radiosondes still show ~100 % duct- 
ing (i.e., the MRG exceeds critical refraction as indicated by 
the vertical dashed-line). The 100-m smoothing shows the 
peak PDF at —420 (N-unitkm _1 ), a very steep gradient that 
is rarely produced in the ECMWF analysis. At the 500-m 
smoothing, most radiosonde profiles are deprived of ducting 
layer. 

On the other hand, the near-coincident ECMWF analy- 
sis (with ~200-m vertical sampling) shows a much broader 
range of MRG variations (from -70 to — 400N-unitkm _1 ) 
that do not match either PDFs of the radiosondes with 200-m 
or 300-m smoothing. The MRG of ECMWF reaches around 
—200 (N-unitkm _1 ) at peak frequency, which is close to 
the radiosondes with 300-m smoothing. But it also shows 
much higher frequency of stronger gradient (<— 240N- 
unitkm -1 ) as compared with the 300-m smoothed radioson- 
des. The large discrepancy between the near-coincident 
ECMWF analysis and the radiosonde with various vertical- 
smoothing levels could not be simply explained by the lim- 
ited vertical resolution in the ECMWF model. Instead, it 
could imply some deficiency in the model physics to repro- 
duce the correct distribution of the MRG at the ABL top. 

Similarly, the PDF of the MRG for COSMIC soundings 
over region (20-25° S, 70-86° W) is plotted against the near- 
coincident ECMWF analysis (Fig. 5b). The ECMWF analy- 
sis shows a similar PDF in Fig. 5a, but with a slightly higher 
frequency for lower MRG, which could be related to weaker 
inversions at higher latitudes towards the south. On the other 
hand, the RO measurements significantly underestimate the 
MRG. The MRG tends to saturate at ^ —1 10 (N-unitkm -1 ), 
which is much smaller than the critical refraction (— 157N- 
unitkm -1 ). 


The radiosondes confirm the prevalence of ducting con- 
ditions at the ABL top over the VOCALS region. In the 
presence of ducting, the standard RO refractivity retrieval be- 
comes a non-unique inversion problem (Sokolovskiy, 2003; 
Xie et aL, 2006). With the assumption of 1-D atmosphere 
with horizontally extended ducting layer at the ABL top, the 
standard RO refractivity retrieval will have a MRG value ex- 
actly equal to the critical refraction, no matter how much the 
MRG exceeds the critical refraction in the real atmosphere. 

The delta function like PDF in the MRG of COSMIC 
(Fig. 5b) implies the dominant presence of ducting, which 
leads to the convergence of the MRG (saturated) to the crit- 
ical refraction threshold in RO soundings. However, the 
weaker refractivity gradients (e.g., — 1 10N-unitkm _1 ) than 
the critical refraction (— 157N-unitkm _1 ) in COSMIC RO 
soundings is partially due to the vertical smoothing (^200- 
m) applied in the standard retrieval process. Moreover, the 
horizontal inhomogeneity in the atmosphere can also intro- 
duce extra smoothing effects (e.g., Ao, 2007) but requires 
further investigation. Overall, the systemically smaller MRG 
in COSMIC soundings should be primarily caused by duct- 
ing but not by the lack of vertical resolution. 

4.3 Refractivity bias (/V-bias) due to ducting in RO 
soundings 

The VOCALS radiosondes and ECMWF analysis confirm a 
wide-spread ducting layer at the top of the ABL over the SE 
Pacific (Fig. 5). The ducting layer has a profound impact 
on the propagation of radio waves and could result in sys- 
tematically negative biases (or called A/ -bias) in RO refrac- 
tivity retrievals (Sokolovskiy, 2003; Ao et aL, 2003; Xie et 
aL, 2006), which result from the non-unique Abel inversion 
problem with details described in Xie et al. (2006). Xie et 
al. (2010a) reveal large A -biases in GPS RO soundings over 
subtropical eastern oceans across the globe, including the SE 


www.atmos-chem-phys.net/12/903/2012/ 


Atmos. Chem. Phys., 12, 903-918, 2012 


910 


F. Xie et al.: Advances and limitations of ABL observations from GPS occultation 


Pacific. The regions with significant N -biases are highly cor- 
related with the area where elevated ducting prevails (von 
Engeln and Teixeira, 2004; Lopez, 2009). 

Here we investigate the PDF of the maximum TV -bias 
caused by the ducting through simple geometric-optics 
(GO) simulations. Based on the VOCALS radiosonde and 
ECMWF refractivity profiles, we conducted a simple two- 
step GO simulation: (1) a forward operator to simulate the 
GPS RO bending angle measurement by integrating an in- 
put refractivity profile, which represents a one-dimensional 
atmosphere; and (2) an inverse operator to simulate the GPS 
RO refractivity retrieval process by integrating the bending 
angle profile (from the first step) through the standard Abel 
integration. The two GO operators that are specifically built 
for simulating ducting cases are fully described in Xie et al. 
(2006). In the absence of ducting, the Abel-retrieved refrac- 
tivity profile will be identical to the input refractivity profile 
and introduce no bias. However, in the presence of ducting, 
the Abel-retrieved refractivity profile becomes negatively bi- 
ased inside the ABL. The solely ducting-induced V-bias (in 
percentile) is defined as (Retrieval - ACinput)//V_i npu t x 100. 
The maximum A/-bias of each simulated profile is generally 
near the top of the ducting layer, i.e. the ABL top. Such 
simulations are applied on all VOCALS radiosondes and the 
near-coincident ECMWF profiles. We then bin the V -biases 
from —20 % to 0 with a 2 % interval. The PDF of the V-bias 
is then normalized by the bin size and the total number of 
samples. 

As shown in Xie (2006), the maximum V-bias introduced 
by the ducting is linearly proportional to the refractivity gra- 
dient at the ABL top and thus is vertical resolution depen- 
dent. In the above GO simulations, we compute the V-bias 
from the radiosonde data with various vertical-smoothing in- 
tervals. Figure 6 shows the PDF of V-bias for radiosonde 
and ECMWF data. Consistent with the results in Fig. 5a, 
the radiosonde profiles with the 100-m and 200-m smooth- 
ing preserve the ducting and the peak frequency of maximum 
V-bias is as large as —12% and —8%, respectively. With 
further smoothing, the V-bias reduces due to a weaker MRG 
in the ducting layer. At the 500-m smoothing, most profiles 
show zero bias, since the MRG rarely exceeds critical refrac- 
tion at this level of smoothing. On the other hand, the PDF of 
the ducting-induced V-bias in the ECMWF is most consis- 
tent with the radiosonde result at the 300-m smoothing, but 
shows a broader spread. 

5 Seasonal mean ABL height climatology and the 
transect along —20° S 

It has been shown above that COSMIC RO is capable of de- 
tecting the ABL top height from the sharp refractivity gradi- 
ent despite the V-bias of RO profiles inside the ABL due to 
ducting. We further generate a seasonal mean ABL height 
climatology using the aforementioned MRG method for 



Fig. 6. PDF of the maximum A -bias due to ducting for radioson- 
des (190 profiles, with different vertical smoothing) and the near- 
coincident ECMWF analysis over the region (15-25° S, 70-86° W). 
The normalized PDF is multiplied by 2 to scale the maximum PDF 
close to 100. 

COSMIC RO soundings over the SE Pacific (0-50° S, 70- 
140° W) during the boreal fall season (September-November 
in 2007, 2008 and 2009). The same method is also ap- 
plied to ECMWF analysis (down-sampled to 1° latitude x 1° 
longitude) profiles in the region. The ABL heights and 
the MRG along the —20° S transect for the VOCALS ra- 
diosondes (200-m vertical smoothing), COSMIC (centered 
at —22.5° S with 5 -deg gridded-average) and ECMWF anal- 
ysis (at 22° S grid) are also analyzed. 

5.1 ABL height climatology 

Deepening of the ABL westward off the coast is clearly ev- 
ident in both COSMIC RO and ECMWF data (Fig. 7a, b). 
COSMIC RO shows an overall deeper ABL and a signifi- 
cant different pattern as compared to the ECMWF. A dipole 
feature in ABL heights is seen in COSMIC RO, i.e., a mini- 
mum ABL height (— 1.2 km) centered at (20° S, 75° W) just 
off the coast of Peru and Chile border; and a maximum 
ABL height (-2.3 km) centered at (15° S, 135° W). The ABL 
heights gradually increase westward (to warmer sea-surface- 
temperature, SST) from a near-coast minimum at — 75° W 
to a maximum at — 135° W within — 8-25° S. A “tongue” 
of low ABL heights extends northwestward from the lo- 
cal minimum and introduces meridional variations in ABL 
heights near the coast. Offshore at higher latitudes (e.g., be- 
yond —25° S and westward of —90° W), the ABL heights 
become rather homogeneous in zonal direction while present 
large meridional gradient (e.g., shallower at higher latitudes), 
which is consistent with the SST pattern (no shown). 

The ECMWF analysis also captures the general feature of 
the deepening ABL westward but shows systematically shal- 
lower ABL (Fig. 7b). A similar dipole feature is also present 
but at different locations, i.e., a minimum (—0.7 km) centered 
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Fig. 7. Mean ABL height climatology (September-November in 2007, 2008 and 2009) for (a) COSMIC RO and (b) ECMWF analysis. 
The standard deviation of the ABL heights is shown in (c) and (d), respectively. The ABL heights and the MRG (mean and standard 
deviation) along the transects at: ~22.5° S (COSMIC), 22° S (ECMWF) and ~20° S (radiosonde) are illustrated in (e) and (f) respectively. 
The horizontal black dashed line in (f) denotes the critical refraction gradient. 


at the Chilean coastal region (^30° S, 72° W) and a maxi- 
mum (~1.7 km) centered at (15° S, 1 15° W). It is also inter- 
esting to note that, ECMWF analysis shows a small merid- 
ional variation in ABL heights between ~8° S and 30° S near 
the coast (70-100° W), but COSMIC RO shows larger merid- 
ional gradients. 

A similar pattern of the ABL height variations (mean 
standard deviation) with comparable magnitude is shown in 
Fig. 7c and d. The ABL height variations tend to be smaller 
over the stratocumulus region near the coast with a local 
minimum of -300-m at (18° S, 85° W) in COSMIC and a 
smaller amplitude (~200-m) slightly northward at (15° S, 
90° W) in ECMWF analysis. The variations increase in both 
westward and meridional directions. A local maximum in 


ECMWF (~0.8km) is clearly seen at ^130° W within 10- 
25° S, where the trade cumulus are located. COSMIC shows 
a slightly smaller maximum at the similar location. Note that 
the trade wind inversion could become decoupled from an 
underneath well-mixed surface layer (e.g., Lock et al., 2000). 
Such a decoupling process could result in double or even 
multiple inversion layers besides the trade wind inversion 
as observed from Hawaii radiosonde soundings (Cao et al., 
2007). These multiple sharp structures result in two-mode 
or even multiple-mode distribution in the refractivity gradi- 
ent PDFs, which could introduce larger variations in ABL 
heights from the MRG method. The large variances in ABL 
heights are also seen over the weak subsidence regions, i.e., 
near the equator (0-5° S) and at higher latitudes near the 
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storm tracks (40-50° S). The weak subsidence tends to pro- 
duce less well-defined ABL and thus results in a broader dis- 
tribution in refractivity gradient PDFs, which could lead to 
noisier ABL heights derived from the MRG method. 

The MRG method performs the best for ABL height de- 
tection when the PDF of the refractivity gradients for each 
profile has one dominant mode, which is generally the case 
over the stratocumulus regions. A separate study on the re- 
fractivity gradient PDFs from COSMIC and ECMWF over 
the trade wind, equatorial (e.g. Intertropical Convergence 
Zone, ITCZ) and high-latitude regions will improve our un- 
derstanding of these more complicated ABL structures and 
help interpret the difference in the ABL height climatology 
between COSMIC and ECMWF. A better and more robust 
ABL height detection algorithm could also be developed and 
benefit from such study. 

5.2 The ABL structure along the —20° S transect 

The bottom two panels in Fig. 7 illustrate ABL height and 
the maximum refractivity gradient along the —20° S tran- 
sect. In COSMIC observation, ABL heights increase steadily 
along the —22.5° S (averaged over 20-25° S with a 5-degree 
longitude bin) transect from off the coast westward to a 
maximum at — 130°W. A relatively steep slope of —30 m 
(100 km) -1 westward is seen from —74° W to — 90°W, 
and the slope decreases significantly further westward to 
be -10 m (100 km) -1 between -100° W and -120° W 
(Fig. 7e). The magnitudes of ABL height variations also 
gradually increases from a near-coast minimum (— 300-m) 
over the stratocumulus regions to a maximum (— 800-m) over 
the trade cumulus regions. 

Similarly, the ECMWF analysis also clearly demonstrates 
the ABL deepening westward along the 22° S transect from 
a near-coast minimum to a local maximum at — 1 05° W 
(Fig. 7e). An identical slope of —30 m (100 km) -1 between 
75° W and 90° W is also seen in the ECMWF analysis, and 
the slope decreases further westward of 90° W before reach- 
ing the local maximum. Contrary to the persistent deepening 
ABL offshore in COSMIC observations, the ECMWF ABL 
heights decrease further westward of — 105° W. Nonetheless, 
the variation of the ABL heights along the transect shows 
slightly smaller amplitude near the coast (— 200-m) but very 
similar amplitude further west in comparison with COSMIC 
observations. 

The ABL heights of individual VOCALS radiosondes 
(15-22° S, 70-86° W) are also shown in Fig. 7 e (red dots). 
The mean and standard deviation of the ABL heights at 
three 5-degree longitude bins (centered at 74° W, 79° W and 
84° W) are also plotted, which matches well with COSMIC 
observations at — 79°W, but slightly higher by —140 m at 
74° W and 84° W. The VOCALS radiosondes also show the 
deepening of ABL from —1.3 km (median value) at 74- 
80° W to —1.6 km at —84° W along the 20° S transect. Such 
radiosonde ABL heights derived from MRG method are con- 


sistent with the ABL heights derived from the radiosonde in- 
version base from the same region (e.g., Fig. 4a of Rahn and 
Garreaud, 2010 and Fig. 10 of Wyant et al., 2010) as well as 
the cloud top height measurements from the airplane cloud 
radar during the VOCALS campaign (e.g., Fig. 2 in Brether- 
ton et al., 2010). 

Figure 7f shows the maximum refractivity gradient at the 
ABL top along the —20° S transect. Note that the MRG is 
vertical resolution dependent. To match the apparent vertical 
sampling of —200-m in ECMWF analysis within —1-2 km, 
the VOCALS radiosonde profiles are vertically smoothed 
with a 200-m window. ECMWF shows a maximum MRG 
near coast (—270 N-unitkm -1 ) and drastically decreasing 
westward to around —110 (N-unitkm -1 ) at — 140° W. The 
variations of MRG also decrease significantly further off- 
shore. The radiosonde soundings also show very similar 
decrease in MRGoffshore but with slightly larger gradients. 
It is important to note that ECMWF analysis shows dom- 
inant ducting presence (MRG < — 157 N-unitkm -1 ) from 
near coast to — 100° W. VOCALS radiosondes also confirm 
the prevailing sharp refractivity gradient. On the other hand, 
the COSMIC RO shows systematically smaller MRG overall, 
with a seasonal mean maximum (-—85 N-unitkm -1 ) near 
the coast and slightly decreasing to -—70 (N-unitkm -1 ) at 
135° W. Again, the MRG bias in COMIC RO can be largely 
attributed to the refractivity retrieval bias due to ducting as 
elaborated in Sect. 4.3. 

5.3 COSMIC RO sampling 

The six-satellite COSMIC constellation is in six inclined or- 
bits with 30° orbit plane separation, which allow full 24-h 
local time sampling within one-month at low and mid lati- 
tudes (e.g., Pirscher et al., 2010; Xie et al., 2010b). Such 
diurnal sampling capability is a very appealing feature for 
compiling the climatology. However, the constellation con- 
figuration also leads to meridional sampling differences, i.e., 
denser sampling per unit area at higher latitudes. It is im- 
portant to understand the sampling errors in the ABL height 
climatology due to the inhomogeneous sampling of COS- 
MIC RO in the meridional direction. As shown in Fig. 8a 
(shaded), during the total of 9-month study period, the sam- 
pling is rather homogeneous in the zonal direction. The num- 
ber of sampling within each 500kmx500km area is —160 
profiles at 40° S and gradually reduces to a minimum of —50 
profiles at 10° S equatorward. Moreover, not all the RO pro- 
files penetrate below 500-ma.m.s.l. The percentage of RO 
penetrating into the lowest 500-m tends to be lower in the 
tropics than at higher latitudes (contours in Fig. 8a). For ex- 
ample, only about 20-30 % RO profiles penetrate deep into 
the lowest 500-m a.m.s.l. near the equator (0-10° S), whereas 
the penetration rate increases at higher latitude and reaches 
over 60 % at —50° S. Note that the penetration rate is rather 
homogeneous in the zonal direction except a slight reduction 
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Fig. 8. (a) Number of COSMIC RO soundings per area of 500kmx500km (shaded); the over-plotted contours denote the percentage of 
soundings that penetrate below 500-m above the m.s.l.; (b) errors in ABL heights caused by the limited COSMIC RO sampling derived from 
ECMWF analysis (contours, with 0. 1 km interval). 


near the coast, where the strongest ducting or temperature 
inversions are located (i.e., within 15-40° S, 70-85° W). 

To estimate potential sampling errors with COSMIC RO, 
we sample the ECMWF analysis at the RO sounding loca- 
tions for the study period and create a sub-sampled ABL 
height climatology (not shown). The difference between the 
full and sub-sampled ABL height climatology reveals the 
sampling errors (Fig. 8b). The sampling errors are generally 
less than 100-m and only become larger at two locations, i.e., 
negative errors in ABL heights up to — 400-m close to (10° S, 
140° W), and smaller positive errors of ~200-m near (35° S, 
130° W). Overall, the limited COSMIC RO sampling results 
in generally small errors in ABL heights and will not affect 
the general morphology of the ABL height climatology in 
Fig. 7a. 

6 Discussions 

Over the SE Pacific, the cool SST along with the large-scale 
subsidence in the free troposphere lead to a well-defined 
ABL with a strong temperature inversion and a sharp mois- 
ture drop across the ABL top. Such a thin transition layer in- 
troduces a very sharp refractivity gradient and thus results in 
a large bending angle in RO signals, which can be precisely 
detected by the high-resolution GPS RO soundings. How- 
ever, there are some limitations in GPS RO for ABL studies, 
which are discussed in the following. 

6.1 Ducting effect on GPS RO 

The sharp refractivity gradient at the ABL top over the VO- 
CALS region often exceeds the critical refraction threshold 
(i.e., —157 N-unitkm -1 ) and results in a ducting condition 
(also referred to as super-refraction). The refractivity gra- 
dient is dominated by the moisture gradient in the tropics 
and subtropics (von Engeln and Teixeira, 2004; Ao, 2007). 


At ~lkm altitude, a mere — 2.5x 10 -2 gkg -1 m _1 in spe- 
cific humidity will lead to a ducting condition. Nearly all 
the VOCALS radiosonde profiles (between 15° S and 25° S) 
meet the ducting condition even after the 100-m smoothing. 
The mean specific humidity gradient at the ABL top of these 
profiles is as large as —5.7 x 10 -2 g kg -1 m -1 with a standard 
deviation of 2.8 x 10 -2 gkg -1 m -1 . 

A horizontally extended ducting layer results in a sharp 
increase in RO bending at the ABL top, which becomes in- 
finitely large (a singularity) in geometric-optics simulation 
(Sokolovskiy, 2003; Xie et al., 2006). By applying the stan- 
dard Abel retrieval (i.e., integrating through the infinite bend- 
ing), the refractivity retrieval based on the spherically sym- 
metric atmosphere assumption produces a maximum refrac- 
tivity gradient of — 157 (N-unitkm -1 ) at the top of the duct- 
ing layer, where the largest bending occurs. The Abel inver- 
sion used for retrieving refractivity becomes a non-unique 
inversion problem and results in negative refractivity biases 
below the ducting layer (Xie et al., 2006), as seen in Figs. 2c 
and 6. The delta-function like PDF of the RO MRG implies 
the prevalence of ducting. The systematically smaller MRG 
in RO refractivity retrieval is primarily caused by the pres- 
ence of a ducting layer at the ABL top. 

6.2 Impact of horizontal inhomogeneous atmospheric 
structure on GPS RO 

Over the VOCALS region, the elevated ducting condition is 
dominant. The idealized 1-D simulation on a horizontally ex- 
tended ducting layer indicates that the maximum gradients in 
RO refractivity retrieval should be close to the critical refrac- 
tion threshold, i.e., —157 N-unitkm _1 (e.g., Sokolovskiy, 
2003; Xie et al., 2006). The maximum bending angle from 
current COSMIC RO measurements rarely exceeds 0.05 ra- 
dians, and so the refractivity gradient rarely exceeds —110 
(N-unitkm _1 ) over the VOCALS region (Fig. 5b). Hor- 
izontal averaging of vertical ABL structures by RO can 
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effectively reduce the observed refractivity gradient (e.g., 
Ao, 2007). In other words, the radiosonde represents a lo- 
calized ABL profile, and the ABL height variability may not 
be uniform over a horizontal span of ^200 km close to the 
RO tangent point, where most GPS RO bending along the ray 
path is accumulated (Kursinski et ah, 1997). Figure 7 clearly 
shows the deepening of the ABL height in both ECMWF 
and COSMIC RO. The slope in the ABL heights (^30 m per 
100 km at ^20° S) indicates the existence of large-scale hor- 
izontal inhomogeneity in the atmosphere. Such smoothing 
effect would need to be further investigated and quantified. 
Also the 200-m vertical smoothing applied on RO bending 
profiles further degrades the vertical resolution when sharp 
bending angles are present. A more sophisticated filter needs 
to be developed to preserve the fine structure of the sharp 
transition in high-resolution bending, which will lead to a 
larger bending and a sharper refractivity gradient (closer to 
critical refraction) at the ABL top. 

6.3 Penetration issue and sampling errors 

COSMIC RO receivers adopt the open-loop (OL) tracking 
technique (Sokolovskiy 2001; Ao et al., 2009), which signif- 
icantly improves the percentage of RO signals penetratating 
deep into the lower troposphere. The OL tracking eliminates 
several problems in previous RO missions (e.g., GPS/MET, 
CHAMP and GRACE) that use the phase-locked loop (PLL) 
tracking technique, i.e., early termination of RO signal track- 
ing and systematic tracking errors caused by the complex 
signal dynamic in the moist lower troposphere (Ao et al., 
2003; Beyerle et al., 2003; Anthes et al., 2008). However, 
the selective lower penetration rate at low latitudes (Fig. 8a) 
in COSMIC sounding requires further investigation. The im- 
proved penetration rate, especially over low latitudes could 
significantly improve the RO sampling and benefit the sci- 
ence application of GPS RO. 

6.4 COSMIC RO ABL height climatology and its 
implication to model analysis 

The difference in the monthly mean ABL heights between 
COSMIC RO and ECMWF analyses (Fig. 7) is substantial 
and cannot be ignored. Such discrepancy is not caused by the 
different vertical resolution and the sampling errors in COS- 
MIC RO and could imply the deficiency in model physics. 
The misrepresentation of MRG and ABL height (Fig. 4) 
imply incorrect strength and altitude of the inversion layer, 
which could have profound impact on the entrainment rate 
and thus cloud evolution simulated in the ECMWF model. 
The climate models with much coarser vertical and horizon- 
tal resolution will have even bigger challenges to correctly 
represent the two critical ABL parameters. 

Due to the lack of ABL observations, the state-of-the-art 
ABL parameterizations are generally derived and validated 
by some field campaign observations. The limited spatial and 


temporal coverage of the measurements could only represent 
the specific local physical processes under certain synoptic 
background forcing conditions. The synoptic and meso-scale 
processes are known to have significant impact on variations 
of the ABL height or vertical structures (e.g., Toniazzo et 
al., 201 1). Therefore the fine-tuned parameterization scheme 
developed from such a dataset could potentially be biased 
to certain atmospheric conditions and limit its general ap- 
plication globally. The global observation of ABL heights 
from COSMIC RO provides an extremely valuable dataset to 
evaluate various ABL parameterization schemes through the 
diagnostic analysis. 

Soon after its launch in 2006, COSMIC RO soundings 
have been operationally assimilated into the leading global 
weather forecasting centers, such as ECMWF and NCER 
Significant positive impacts of RO observations on the upper 
troposphere and lower stratosphere have been demonstrated 
(Cucurull et al., 2006; Healy and Thepaut 2006; Cucurull et 
al., 2007). However, the much-anticipated positive impact of 
RO soundings in the lower troposphere, especially over the 
mid and low latitudes has yet to be demonstrated. 

NCEP operationally assimilates COSMIC refractivity pro- 
files. However, the systematic A -bias in RO refractivity over 
the subtropics eastern oceans (Xie et al., 2010a) along with 
the complicated error characteristics in the lower troposphere 
results in very limited impact of RO soundings on model 
analysis. 

ECMWF, on the other hand, assimilates RO bending an- 
gles, which should not be biased due to ducting. However, in 
the presence of ducting, RO bending angle loses information 
inside the ducting layer and a layer right below, and becomes 
a singularity (e.g., an infinite bending angle) at the top of the 
ducting layer (Sokolovskiy, 2003; Xie et al., 2006). More- 
over, directly assimilating RO bending angle in the presence 
of ducting will not avert the non-unique inversion problem, 
i.e., a single bending angle corresponds to a continuum of 
refractivity profiles (Xie et al., 2006). For example, given 
the constraint of the observed RO bending angle profile, the 
forward calculation (e.g., computing bending from model re- 
fractivity) will not tell the difference among the continuum of 
model refractivity (or temperature, humidity) profiles with or 
without the bias. Therefore, to maximize the benefit of the 
fine vertical resolution RO observations in the lower tropo- 
sphere for global models, the data assimilation system needs 
to develop the capability to deal with the singularity problem 
due to the ducting. 

7 Conclusions and future work 

In this paper, we analyzed the ABL structures in COSMIC 
GPS RO, radiosondes and ECMWF analysis data over the 
SE Pacific where the VOCALS campaign took place. The 
VOCALS radiosonde observations show a well-defined ABL 
top with a strong temperature inversion and a sharp moisture 
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gradient, which produces a sharp refractivity gradient and a 
large bending angle in the GPS RO signals. The refractivity 
gradients over the region are often so large that they exceed 
critical refraction (e.g., —157 N-unitkm _1 ) and become the 
ducting (or super-refraction) condition. This situation can re- 
sult in a systematic negative refractivity bias (TV -bias) in the 
standard RO retrieval inside the ABL (Sokolovskiy, 2003; 
Xie et al., 2006). It is noteworthy that the A/-bias does not af- 
fect the ABL height detected with the MRG method. Instead, 
the very sharp refractivity gradient caused by ducting allows 
the best detection of ABL height from GPS RO through the 
simple MRG method. 

The comparison of RO ABL height derived from the MRG 
method with radiosondes is encouraging, however, ECMWF 
shows a systematically low bias in the ABL height and the 
bias becomes larger in the deeper ABL (Fig. 3b). The VO- 
CALS radiosondes exhibit a wide range of ABL height vari- 
ations (1^2 km) in the PDF, with the peak occurrence fre- 
quency at ~1.3 km (Fig. 4b). The vertical smoothing results 
in a slightly lower ABL height, but does not change the over- 
all shape of the PDF. The near-coincident COSMIC RO cap- 
tures the wide range of ABL height variations (Fig. 4a, c). 
ECMWF analysis, on the other hand, tends to produce sys- 
tematically shallower ABL as compared with both the near- 
coincident radiosondes (Fig. 4b) and COSMIC RO sound- 
ings (Fig. 4c). The discrepancy between the ECMWF and ra- 
diosondes is not due to the limited vertical resolution (^200- 
m) in the model. 

The PDF of the radiosonde MRG depends strongly on the 
vertical resolution of the refractivity profile. The higher or- 
der smoothing results in a smaller MRG value. Nearly all 
the VOCALS radiosonde data have a ducting layer even af- 
ter 200-m smoothing. The near-coincident ECMWF analy- 
sis underestimates the high MRG and overestimates the low 
MRG as compared to radiosondes after 200-m smoothing. 
On the other hand, COSMIC RO shows a systematic bias in 
the MRG and the RO refractivity inside the ABL, which are 
primarily caused by the non-unique inversion problem in the 
presence of ducting. The ABL height derived from COS- 
MIC RO is generally not affected by the ducting. However, 
it would require specific attention for directly using such bi- 
ased RO refractivity profiles for scientific studies. Simulation 
study based on radiosonde soundings reveals strong depen- 
dence of ducting-induced A-bias on vertical resolutions, i.e., 
higher order smoothing leads to a smaller A -bias (Fig. 6). 
The 100-m smoothing shows the peak frequency of A-bias 
~ — 12 %, which reduces to ~ — 8 % with 200-m smoothing 
and decreases further with additional smoothing. ECMWF 
analysis shows the PDF of A-bias comparable to the ra- 
diosonde soundings with ~300-m smoothing, but underes- 
timates the ducting occurrence frequency. 

Both COSMIC RO and ECMWF analysis clearly show 
the deepening of the ABL heights westward from the shal- 
low stratocumuli near the coast of Peru and Chile to the 
much deeper trade wind cumuli. However, the difference 


in ABL height between ECMWF analysis and COSMIC RO 
is significant and cannot be ignored. COSMIC RO shows 
a deeper ABL with a minimum ABL height (~1.2 km) cen- 
tered at (20° S, 75° W), which is located ~10° north of what 
is seen in ECMWF analysis. The maximum ABL height 
(~2.3 km) from COSMIC RO is centered at 20° west of 
ECMWF (Fig. 7). At low latitudes, despite the decreasing 
number of COSMIC RO soundings and the lower percentage 
of soundings that penetrate into the lowest 500-ma.m.s.l., 
there are small sampling errors in the mean ABL height cli- 
matology. 

The ABL height observed from COSMIC RO provides 
an extremely valuable and independent dataset for diagnos- 
tic analysis of various ABL parameterization schemes in 
weather and climate models. The biased RO refractivity due 
to the ducting over the eastern oceans and the trade wind re- 
gions (Xie et al., 2010a) would need to be corrected (e.g., 
Xie et al., 2006) before scientific application. Improving the 
percentage of RO profiles penetrating into the ABL and a bet- 
ter understanding of horizontal averaging effect in RO would 
further enhance the scientific values of GPS RO soundings 
for global ABL studies. 
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